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We discuss the cosmic ray spectra in annihilating/decaying Nambu-Goldstone dark matter models. 
The recent observed positron/electron excesses at PAMELA and Fermi experiments are well fitted 
- - . by the dark matter with a mass of 3 TeV for the annihilating model, while with a mass of 6 TeV 

for the decaying model. We also show that the Nambu-Goldstone dark matter models predict a 
I distinctive gamma-ray spectrum in a certain parameter space. 

(N : 

bJO' I. INTRODUCTION 

^ : 

! Recent observations of the PAMELA [l| , ATIC [1| , PPB-BETS , and Fermi Q experiments strongly suggest the 
' existence of a new source of positron/electron fluxes. The most interesting candidate of the new source is the 
fS| ' dark matter with a mass in the TeV range annihilating or decaying into the visible particles which result in the 
positrons/electrons [5] ■ Especially, the Fermi experiment has released data on the electron/positron spectrum from 
20 GeV up to 1 TeV where the spectrum falls as E~^-^. As reported in Refs. [1,0, Hi the data can be well fitted by 
^. the dark matter which mainly annihilates/decays into a pair of light scalars each of which subsequently decays into 
I ' a pair of electrons or muons.^ 

From the theoretical point of view, it is always motivated to relate the identity of the dark matter with the new 
physics which is anticipated from other motivations [ll[ . Among them, one of the most motivated new physics is the 
supersymmetric standard model (SSM) which is expected as a solution to the hierarchy problem. In the SSM, the 
stable hghtest supersymmetric particle (LSP) is a candidate of the dark matter. The LSP interpretation with a mass 
in the TeV range [12, JSJ, however, implies that the masses of the other supersymmetric particles are much heavier 
than a TeV, which diminishes the significance of the SSM as a successful solution to the hierarchy problem. Besides, 
such a heavy LSP interpretation can be falsified rather easily, once the supersymmetric particles are discovered in 
hundreds GeV range at the coming LHC experiments. 
fT^ ' The attempt to relate the dark matter to the supersymmetric models, however, should not be necessarily confined 
, to the LSP dark matter scenarios. In fact, the SSM always requires other new physics, the supersymmetry breaking 
■ sector which may include a stable particle as a candidate of the dark matter. The idea of the Nambu-Goldstone dark 
^\ matter developed in Ref. is one of the realization of the dark matter in a supersymmetry breaking sector. There, 
. the dark matter is interpreted as a pseudo-Nambu-Goldstone boson in a supersymmetry breaking sector so that the 
^ ' mass of the dark matter is in the TeV range out of the dynamics of supersymmetry breaking at tens to hundreds TeV 
range. 

In an explicit example of the Nambu-Goldstone dark matter scenario given in Ref. we considered the model 
where the dark matter annihilates into a pair of the light pseudo scalars, the R-axions, via a narrow resonance, the 
flaton, which leads to the right amount of the dark matter. Furthermore, for the R-axion mainly decaying into a pair 
of electrons or muons, the final state of the dark matter is the four electrons or muons, which are favored to explain 
the positron/electron spectrum observed by Fermi experiment. 

The most prominent difference of the Nambu-Goldstone dark matter model from the other dark matter models 
which explain the observed positron/electron excesses is that the Nambu-Goldstone dark matter model is strongly 
interrelated to the physics of the SSM. That is, in the Nambu-Goldstone dark matter scenario, it is difficult for the 
supersymmetry breaking to be much higher than tens to hundreds TeV. This restriction suggests that the supersym- 
metry breaking effects should be mediated to the SSM sector at the low energy scale, i.e., the model requires the gauge 
mediation mechanism ^1^. Thus, the dark matter interactions with the SSM particles are determined along with the 
gauge mediation effects. With the interrelation to the SSM physics, the Nambu-Goldstone dark matter model has a 
distinctive prediction on such as a gamma-ray spectrum and an antiproton flux in cosmic ray. 
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^ The earlier works on the cosmic ray spectra before the Fermi data in the presence of the light scalars decaying into the light lepton pairs 
can be found in Ref. [9l[lo|. 
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In this paper, we discuss the cosmic ray spectra in the Nambu-Goldstone dark matter scenario based on the 
exphcit model in Ref. [ij]. As we will show, the Nambu-Goldstone dark matter model fits to the recently observed 
positron/electron excesses in both the annihilating and the decaying dark matter scenarios. Furthermore, we show 
that the model gives a distinctive prediction on a gamma-ray spectrum for a certain parameter space, which comes 
from the finite annihilation/decay rates of the dark matter into a pair of gluinos. 

The organization of the paper is as follows. In scctionlTTl we summarize the explicit example of the NDGM model 
in Ref. [l^. There, we also derive conditions so that the dark matter density is consistent with the observed dark 
matter density. In section llTTl we show how well the model fits the observed fluxes in both the annihilating/decaying 
dark matter scenarios. We also demonstrate how the modes into the SSM particles affect the gamma-ray spectrum. 



II. NAMBU-GOLDSTONE DARK MATTER 



In this section, we review a model of the Nambu-Goldstone dark matter model which is based on a vector-like 
SUSY breaking model in Ref. In the model, the dark matters are interpreted as pseudo Nambu-Goldstone bosons 
which result from spontaneous breaking of the approximate global symmetry in the vector-like SUSY breaking model. 



A. Nambu-Goldstone Dark Matter, Flaton, and R-axion 



The key ingredients of the model are the light flaton which corresponds to the so-called pseudo moduli of the SUSY 
breaking model, and the R-axion which is a pseudo Nambu-Goldstone boson resulting from spontaneous R-symmetry 
breaking. In the foUowings, we overview the relevant properties of those particles. 



1. SUSY breaking sector 

The vector-like SUSY breaking model is based on an SU{2) gauge theory with four fundamental representation 
fields Qi{i — 1, • • • ,4) and six singlet fields Sij — ~Sji — 1, • • • ,4). In this model, the SUSY is dynamically 
broken when the Q's and S"s couple in the supcrpotential, 

W^K,S,,Q,Qj,{i<j), (1) 

where Ay denote coupling constants. The maximal global symmetry this model may have is SU{4) ~ S0(6) symmetry 
which requires A^ = A. The SUSY is broken as a result of the tension between the F-term conditions of S"s and Q's. 
That is, the F-term conditions of Sij, dW/dSij — XijQiQj — 0, contradict with the quantum modified constraint 
Pf(My ) = A^yjj where denote composite gauge singlets made from QiQj. 

Below the dynamical scale Adyn, the model is described by the fight degrees of freedom, and Sa, {A — 1 — 6), 
with a quantum modified constraint, 

J2 Ml^Aly^. (2) 

A=l-6 

Here, we have assumed that the effective composite operators Ma are canonically normalized. Notice that the above 
quantum modified constraint breaks the global SU{A) ~ 50(6) symmetry into SP(2) ~ 5*0(5) symmetry. Thus, when 
the model possesses the 50(6) symmetry approximately, the Nambu-Goldstone bosons appear as the result of the 
quantum modified constraint. Furthermore, the Nambu-Goldstone bosons have a long lifetime when an appropriate 
subgroup of 50(5) is almost exact. (We will discuss more on the stability and the global symmetry in the next 
section.) In this way, we realize a model with the Nambu-Goldstone dark matter with mass in the TeV range out of 
physics of tens to hundreds TeV range. 

To make the discussion concrete, for a while, let us assume that the SUSY breaking sector possesses an 50(5) C 
50(6) global symmetry and rearrange the tree-level interaction Eq. ([1]) so that the 50(5) symmetry is manifest; 

W = A5o(gQ)o + A' ^ Sa{QQ)a, 

a=l-5 

= AAdy„5oA/o + A' J2 ^dynSaMa . (3) 
a=l-5 
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In the second line, we rewrite the superpotential by using the low energy field M^.^ We further assume A < A' 
and assume that the model possesses the iS'0(6) symmetry in the limit of A ^ A'. Under these assumption, the 
Nambu-Goldstone bosons are charged under the SO{5) symmetry and are stable. 
For later convenience, let us solve the quantum modified constraint explicitly by, 



a=l~5 



and plug it into the effective superpotential in Eq. ^ ; 



Weff ~ AA2y„5o- ^SoMl+ E A'Ad.y„5aMa + 0(M4) . (5) 



2 

a— 1 — 5 a— 1 — 5 



From this expression, we see that the SUSY breaking vacuum is given by, 

Fs„ = AA2y„ , Sa^O, Ma^O. (6) 



2. Mass spectrum of So multiplet 

Around the vacuum in Eq. ([6|), the tree-level potential in the 5*0 direction is fiat, and the masses of the Sq multiplets 
can be significantly lighter than the dynamical scale. ^ As discussed in Ref. [3], this is a quite favorable feature to 
account for the observed dark matter density. The annihilation cross section of the Nambu-Goldstone dark matter 
with a mass in the TeV range is generically too small to explain the observed dark matter density (see Ref. fli*! for 
general discussion). In this model, however, the annihilation cross section is enhanced by a narrow resonance which 
is served by the fiaton, a radial component of Sq. It should be noted that, for that purpose, the R-symmetry must be 
broken so that the dark matters with R-charge zero annihilate through the flaton resonance with R-charge two. 

In Ref. [l^ . we considered the model with spontaneous R-symmetry breaking. The spontaneous R-symmetry 
breaking is caused by effects of higher dimensional operators of So in the Kahler potential. 
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where A's denote the dimensionful parameters and the ellipsis denotes the higher dimensional terms of Sq. The above 
Kahler potential provides an effective description of a quite general class of the models with spontaneous breaking of 
the R-symmctry breaking. Especially, when the above Kahler potential results from radiative corrections from physics 
at the scale Adyn, the dimensionful parameters are expected to be, 

Al 167r2Ai^„' A| IGn^Ai^^' 

where dimensionless coefficients C4,6 are of the order of unity (sec Ref. 14] for an explicit model with R-symmetry 
breaking). From this Kahler potential, the R-symmetry is spontaneously broken by the vacuum expectation value of 
the scalar component of Sq; 

where we have introduced the R-symmetry breaking scale fn ^ O(Adyn). 

Now let us consider the masses of the scalar components of Sq. At this vacuum, the scalar component of Sq is 
decomposed into the fiaton s and the R-axion a by, 

^o = ^(/« + s)e-//-. (10) 



^ In the second line, we neglected order one coefRcients of each term. 

^ The fermion components of So corresponds to the spin one half component of the gravitino. 
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Then, the mass of the flaton is given by, 

where we have used — AA^^^^ and assumed Eq. ([5]) with C4 = cg = 0(1). Therefore, the flaton can be in the TeV 
range for A ~ 1 and C4 ~ 1, which is a crucial property for the flaton to make the narrow resonance appropriate for 
the dark matter annihilation.^ 

The R-axion mass, on the other hand, can be much more suppressed, since it is a pseudo-Goldstone boson of the 
spontaneous R-symmetry breaking. For example, when explicit breaking of the R-symmetry mainly comes from the 
constant term in the superpotential,^ the R-axion acquires a small mass [17j . 

mLon - (12) 

JR 

This mass is in a MeV range for fji ~ \/Fso — 30 TeV and TO3/2 — 1 eV, for example. 

The most interesting feature of the R-axion in the above mass range is that it mainly decays into an electron 
pair [Tsl . Therefore, the dark matter which annihilates (or decays) into a pair of the R-axions ends up with the flnal 
states with four electrons. In section lllll we see that the four electron final states of the dark matter annihilation or 
decay are favorable to explain the positron/electron spectra. Notice that there is a lower limit on the decay constant 
fn of the axion- like particles which mainly decay into pairs of electrons from the beam-dump experiments p^oj ; 

(jji \ 1/2 
— ^ for 1 MeV <ma< 100 MeV. (13) 

IMeV/ a K J 

On the other hand, as we will see in section HTC] the right amount of dark matter requires /r ~ 30 TeV. Thus, the 
R-axion decaying into a electron pair in Nambu-Goldstone dark matter models is consistent with the constraint from 
the beam-dump experiment (see for example Ref. ^2d\ for detailed constraints on the axion-like particles). 

One may also consider explicit R-symmetry breaking by higher dimensional terms which are suppressed by some 
mass scale, M*, such as a quartic term Sq/M^..^ With this breaking term, the R-axion mass in hundreds MeV range 
is realized for ffi ~ 30 TeV and M» ~ 10^®~^^ GeV. In this case, the R-axion mainly decays into a pair of muons. In 
the following analyses, we consider the R-axion of a mass in both a MeV range and hundreds MeV range. 



3. Mass spectrum of Sa and Ala multiplets 

We now consider the spectrum of other light particles, Sa and Ma- From the superpotential Eq. ([5]), the fermion 
components of Sa and Ma obtain masses of O(AAdyn)- Scalar masses, on the other hand, comes from the scalar 
potential. 



V = |A'AdynAfaP + 



1 



AJi-dyn 

A2 



\XSoMa + X'AdynSa\ 



- (A- - A^)A^,„ + -(. + fnY 4, + ^ (A- + A^)A^,„ + -is + y, 



A2 



AA' 



dyn 

AA' 



+ 7^^' Kyn^s + ;^A'-A5y„y^ + -^Ady„(s -I- fn) XmXs + -^Adyn(s -I- Jr) 



2 4 A^ 2 2 2 

+ A Adyn -I- {x^ + Vm) J 



V2 



(14) 



Even for large couplings, A, 04^61 there is a possibility that the model involves a "light" flaton of a mass in the TeV range, although 
those models are incalculable. 

^ In this study, we assume that the messenger sector of the gauge mediation also respects the R-symmetry. Otherwise, the radiative 
correction to the Kahler potential of Sq from the messenger sector gives rise to the dominant contribution to the R-axion mass. The 
R-breaking mass from the Higgs sector, on the other hand, is smaller than the one in Eq. I|12|l . even if the so-called fi-term does not 
respect the R-symmetry. 

® Here, we assume Zgn symmetry instead of the continuous R-symmetry. The main conclusion of this paper is not affected as long as the 
order of the discrete symmetry is high enough. 
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where we have decomposed the scalars by, 

Ma = -^{Xm + iy,n), (15) 

and suppressed the index of 5*0(5). Notice that the R-axion does not show up in the scalar interactions in this basis, 
and it only appears in the derivative couplings. 

From the above potential, we find that the pseudo-Nambu-Goldstone mode resides not in {ym,ys) but in [xm,Xs) 
and the eigenvalues of the squared mass matrix of {xm, Xg) are given by, 

\ (trM^ - V(trAf^)^-4detM^) = ^ . ^ (16) 

ml = \ (tr Af2 + V(trA/2)2_4detM2) ^ \^ (a^^^ + {S,f) , (17) 

trM2 = (2A'2 - A2)A2y,, + A2 (5o)' , (18) 
detM^ = A'2(A'2 - A2)A4y^, (19) 

where the rightmost expressions of m| ^ are valid for A ~ A'. The eigen mode </> corresponds to the pseudo-Nambu- 
Gladstone boson of the approximate 5*0(6) symmetry, which becomes massless in the limit of an exact 5*0(6), i.e. 
A A'. The masses of (y™, ys) are, on the other hand, of O(AAdyn)- 

The R-axion interactions only appear in the kinetic terms. In the current basis, the R-axion interactions come from 
the kinetic terms of Sq and Sa, 

C^\{daf(l + ^\ +^{daf(xl+vl) + ^d^a{xsd^ys-ysd^xs). (20) 

^ V JrJ ^jr jr 

Altogether, in the Nambu-Goldstone dark matter scenario (i.e.. A' — A <C 1), the light particle sector below 
0(10) TeV consists of the dark matter and the flaton in the TeV range, and the gravitino and the R-axion with much 
smaller masses, while the other components of Sa and Ma have masses of the order of the SUSY breaking scale, 
A^/'^Adyn- The most relevant terms for the dark matter annihilation are, then, given by, 

where the first term comes from the scalar potential in Eq. (jl4p . while the second term comes from Eq. (j20p . 



4- Flaton decay 

Now let us discuss the decay of fiaton which is important to estimate the dark matter annihilation via the s-channel 
exchange of the flaton. First, we consider the decay mode into a pair of the R-axions. The relevant interactions for 
the decay come form the first term in Eq. (|2ip , and the decay rate into a pair of the R-axion is given by. 



I ml . / \2 /30TeV^ ^ 



r._M = - 4 X 10- j j , (22) 

where we have neglected the mass of the R-axion and taken the final state velocity to be /?/ = 1. 

Next, we consider the flaton decay into a pair of the dark matter. The relevant interaction term is given in Eq. (j21|) 
and the resultant decay rate is given by, 

where ji^ denotes the size of the velocity of the dark matter. Notice that the value of Y s^^^j (i^ is well-defined even 
in the unphysical region, i.e., 2m^ > mg. As a result, we find that Fg^^^//?^ is suppressed compared with Fg^aa, 
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where we have used ~ and uih ^ m^. 

As we will see in section llll 5[ the flaton also decays into a pair of the SSM particles, which are roughly suppressed 
by the mass ratio squared of the SSM fields and the flaton compared with the R-axion mode. Thus, to determine the 
dark matter density, it is good enough to consider the decay into the R-axion. Since we are mainly interested in the 
R-axion of a mass in a MeV or hundreds MeV range, the R-axion in the final state eventually decays into pairs of 
electrons or muons. Therefore, the dark matter annihilation ends up with the final state with four electrons or four 
muons. 

Put it all together, we obtain the flaton decay width at i?cM — tUs, 

^s{Ecm) = Ts-^aa + Fs— 00 -f • • • . (25) 

where -Ecm > Sm^, and the ellipses denotes the modes into the SSM particles (see section lTlI 5p . In the following 
analysis, we approximate the above decay rate by, 

TsiEcu) ^Ts{ms) C^Ts^aa, (26) 

since the other modes are subdominant for Eqm — ms- 

B. Symmetry and Stability of Dark Matter 

So far, we have assumed that the SO{5) global symmetry out of the maximal SO{6) is exact and the dark matter 
is completely stable. This symmetry has been imposed by hand to make the dark matter stable. Although there is 
nothing wrong with this assumption, it is more attractive if the stability of the dark matter is assured by symmetries 
which are imposed by other reasons than the stability of the dark matter. In this section, we show that the present 
model can have such an accidental symmetry by which the stability of the dark matter is achieved. Once the stability 
is assured by an accidental symmetry, the dark matter decays only through higher dimensional interactions which 
do not respect the accidental symmetry. As we will see in the following sections, the lifetime of the dark matter is 
long enough, and furthermore, the lifetime can be in an appropriate range to explain the observed positron/electron 
excesses. 



1. Stability and accidental symmetry 

In the SUSY breaking model discussed above, we may take a U{1) subgroup of the global S'C/(4) symmetry a gauge 
symmetry. As discussed in Ref. the SUSY breaking model with such a gauge symmetry ensures spontaneous 
R-symmetry breaking in certain parameter space (see also Ref. [2]|).^ 

We assign the U{1) charges to the fundamental fields; Qi{l/2), (52(1/2), Q3(— 1/2) and Q4(— 1/2), which corre- 
sponds to; Mi2(l), M34(— 1), while other mesons are neutral. The charges of the singlets S"s are assigned so that the 
interaction in Eq. ([T]) is consistent with the U{1) symmetry. In this case, the low-energy effective superpotential of 
the gauged lYIT model is given by 

W = XS+M^ + XS-M+ + XaSaMa , (27) 

with the quantum constraint 2M_|_M_ + MaMa — ^\yn- Here, the subscript a runs a = 1 — 4, and M± corresponds to 
Mi2 and M34, respectively. For A < Aa, the quantum constraint is satisfied by 2M+Af_ = A^y^^, and the [/(I) gauge 
symmetry is broken down spontaneously. 

Notably, with just the introduction of the U{1) gauge symmetry, the model now has an accidental symmetry which 
assures the stability of the dark matter under the gauge symmetries and the R-symmetry. Here, let us remind ourselves 
that none of these interactions are imposed to assure the stability of the dark matter. 

Under the above symmetries, the lowest dimensional interactions in the superpotential are those in Eq. ([1]), and the 
second lowest dimensional interactions are suppressed by some high mass scale M^, such as terms in the superpotential, 
SQ'^/M^. Here, we are assuming the R-charge assignment; 5(2) and Q{0) (see discussion at the end of this section). 
As a result, the model possesses an accidental global C/(l) symmetry, at the energy scale much lower than M,, with 



^ The introduction of the U{1) gauge symmetry solves the domain wall problem in the original vector-like SUSY breaking model in 
Ref. [3. 
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the charge assignment; S{2) and Q{—1), which is in terms of the low energy fields; S{2) and M(— 2). This symmetry 
is broken down to a discrete symmetry by an anomaly to the SP{1) dynamics. Altogether, in the low energy 
theory, the model has an accidental Z4 symmetry (which is effectively a Z2 symmetry in terms of the low energy 
effective fields, S's and Af's). 

After spontaneous SUSY breaking, the R-symmetry and the U{1) gauge symmetry are broken spontaneously by 
the vacuum condition 2Af+Af_ = A^^^^ and S+S- 7^ 0. Then, an appropriate combination of the accidental Z4 
symmetry and the U{1) gauge symmetry remains unbroken at the above vacuum. Concretely, a rotation under the 
U{1) gauge symmetry by an angle tt on M± is cancelled by a change of signs of M± under the Z4 symmetry, while 
the dark matters change signs under this unbroken combination. Therefore, there is an accidental Z4 symmetry at 
the vacuum, under which the dark matter changes sign. 

Put it all together, in the SUSY breaking model with the U{1) gauge symmetry, the dark matter stability is 
assured by the accidental Z4 symmetry which results from the gauge symmetries and the R-symmetry. The lowest 
dimensional interactions which break the Z4 symmetry are suppressed by some high mass scale M^, and hence, the 
lifetime of the dark matter decaying via the symmetry breaking interactions is much longer than the age of the 
universe. Furthermore, as we will see below, the lifetime of the dark matter by the decay via the lowest dimensional 
symmetry breaking terms is in an appropriate range to explain the observed positron/electron spectra. 



2. Decaying dark matter 



The lowest dimensional interactions which violate the accidental Z4 symmetry are, for example, given by the 
following terms in Kahler potential, 

K = CoY^ Q^Q3 = C Adyn ^ Ma, (28) 

a=l-4 

which respect gauge symmetries as well as the R-symmetry. Here, we have omitted flavor indices of dimensionless 
coefficients, Cq and c, and rewritten the operators in terms of the low energy fields in the rightmost expression. Notice 
that the above operators have physical effects only through supergravity, and hence, the effects are suppressed by the 
Planck scale Af, ~ AfpL ~ 2.4 x 10^^ GeV. 

In supergravity, the above terms result in mass mixings between the flaton and the Nambu-Goldstone dark matters, 

C / Adyn \ I fli \ f\ \ \2 _ JC 2 



^dyn 



^^"VlV^j ) (^^dyn)' X„S = W X^S, (29) 



and the mixing angles are given by. 



3m 



2 ' 



(30) 



where we have omitted flavor indices of the dark matters again and used ~ 2to0 in the final expression. With 
these mixings, the dark matters decay in a similar way of the flaton discussed in section lll A41 ^ and the main decay 
mode is that into the R-axion pair and the lifetime is given by. 



Here, we have taken ~ A^yn for simplicity. As a result, the lifetime of the dark matter is much longer than the 
age of the universe. 

Before closing this section, let us comment on how accidental the above "accidental" symmetry is. In the above 
discussion, we have assigned R-charges S{2) and Q(0), so that the R-symmetry forbids the Z4 breaking terms sup- 
pressed by a factor of Af~^, while allowing the breaking terms suppressed by Af~^ as given in Eq. ([28|) . With this 
charge assignment, however, linear terms of Sa in a superpotential are still consistent with all the other symmetries 
than the "accidental" Z4, while they break the "accidental" Z4 symmetry. Thus, strictly speaking, the Z4 symmetry 



With the above symmetry breaking term, the dark matter decays into the MSSM particles with the same suppression factor Mpj^. The 
branching ratios of those modes, however, are of 0(m?/A® ) and highly suppressed. 
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is hard broken, and hence, it cannot be an accidental symmetry. To avoid this problem, one may consider another 
charge assignment of the R-symmetry; S{2/5) and (3(4/5). Under this charge assignment, the linear terms of 5"s are 
forbidden and the lowest Z4 breaking terms are given by /M"^, which result in similar mass mixings between the 
dark matter and the flaton given in Eq. (|30p . In this way, we may have a truly accidental Z4 symmetry which is the 
outcome of the gauge symmetries and the R-symmetry.^ 

C. Dark matter density 

In this section, we discuss the parameter space which reproduces the observed dark matter density by the Nambu- 
Goldstone dark matter. 

1. Annihilation cross section via flaton resonance 

As derived in Ref. [3], the annihilation cross section of the dark matter into a pair of the R-axion via the flaton 
resonance is given by, 

327r /3, \ml - m\ ) {El^ - m^)^ + m^^l " 647r 1^ A^^,, + 2/| ) m| (<5 + v^^J^f + 7?' ^ ^ 

where the and are defined at -Ecm > 2m0. In the final expression, we have used the non-relativistic approxi- 
mation, 



^Im = + m\vl^, , /30 = y'l - ^vfillEl^ , (33) 
and introduced parameters b and 7s by, 

ml^^m\{\-b), 7s = rs/ms. (34) 
Notice that the cross section does not depend on the parameter A explicitly. 

2. Required cross section 

In the presence of the narrow resonance, the thermal history of the dark matter density is drastically changed from 
the usual thermal relic density without the resonance |;22i, |23|. As a result, the required annihilation cross section 
to account for the observed dark matter density [23 | is different from the one for the non-resonant annihilation cross 
section, 

{av,^i) - 10-9 GeV. (35) 

Instead, in terms of the annihilation cross section at the zero temperature, the required annihilation cross section to 
obtain the correct abundance is given by (25| , 

{cTv,,i) |t=o - 10-9 GeV-2 X — . (36) 

Here ~ 30 denotes the freeze-out parameter of the usual (non-resonant) thermal freeze-out history, while the 
effective freeze-out parameter Xb is defined by, 

Xb (O-Wrol) |t=0 Jxf 



As an alternative solution, we may forbid the linear terms of S's by assuming a conformal symmetry in the limit of the vanishing 
gravitational interactions. Once the linear terms of S's are forbidden at the tree-level, they are not generated by any radiative corrections. 
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For a well tuned and very narrow flaton, i.e. |<5|,7s 1, we found that the effective freeze-out parameter is fairly 
approximated by, 



(38) 





■ 5 ' 




arctan 




) 


V2 


.7s. 





From this expression, we find the asymptotic behaviors of the effective freeze-out parameter as follows. 



Unphysical pole ((5 > 0), (5 < 7s < 10 



-1 



• Unphysical pole ((5 > 0), 7^ < 5 < 10"^ 

• Physical pole {5 < 0), \S\ < 7s < 10"^ 

• Physical pole (5 < 0), 7s < \S\ < 10"^ 



^ ~ 7r7s . (39) 
a;^^~2J. (40) 
x^i~^7s. (41) 



x^^ ~ 27r — . (42) 

Here, we have used arctan(l/a;) — 7r/2 — arctan(a:) for x > in Eq. ([^0]) . and arctan(l/x) = —tt/2 — arctan(x) for 
X < in Eq. (|40p . Notice that the Breit-Wigner enhancement is realized for the first three cases, that is, 



(fTVrei) \t=o » 10-' GeV-2 . (43) 
while the late time cross section can be smaller than the non-resonant annihilation in the final case, i.e., 

(fTz;,.ei)|T=o« 10-'GeV-2 . (44) 



3. Constraints on SUSY breaking scale 

Now let us compare the cross section of the Nambu-Goldstone dark matter via the narrow flaton resonance in 
Eq. ([5^ and the required cross section in Eq. ([M]) . In the hmit of \6\ <C 7s (see Eqs. ([M)) and (HI])), the required 
annihilation cross section in Eq. (j36p is reduced to, 

1 f2 

(CTWroi) \t=o ~ 10^' GeV-^ X ~ 8 X 10"' GeV"^ x . (45) 

x/7r7s Xf TO^ 

Here, we have used the flaton decay rate in Eq. ([22]) . In the same limit, the predicted cross section in Eq. ((32|) is 



where, again, we have used 7s in Eq. (j22[) . 

First, we consider the case of the unphysical pole, 5 > 0. In this case, the require cross section decreases in S^^ 
for a larger S and in a region of 6 > js, while the predicted cross section decreases in S~'^. Therefore, in order for 
the predicted cross section meets the required value, the predicted cross section in the limit of 5 ^ 7s must be larger 
than the required one, i.e., 

''^^ > 8 X 10-9 GeV-2 X (47) 



Thus, for X/ ~ 30, we find a constraint on //j and Adyn 



/fl,< 50TcVx ( l + :^j . (48) 



Notice that this bound is independent of the mass of the dark matter. 
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10"^ 10"' 10"^ 10"^ 10"^ 10"^ 10"* 10"^ 10"^ 10"^ 



6 6 

FIG. 1: Left) The S dependence of the required annihilation cross section at zero temperature from the observed dark matter 
density in the case of the unphysical pole (red line) for js — A x 10~*. = 3TeV, Adyn = 30TeV and /_r = 30TeV, which 
corresponds to The green line shows the required annihilation cross section in the usual thermal history. The blue line shows 
the predicted annihilation cross section for = 3 TeV, Adyn = 30 TeV and fn — 30 TeV, which corresponds to 7s = 4 x 10~*. 
Right) The required annihilation cross section in the case of the physical pole. In the figures, we have assumed Xf = 30. 



In Fig.m we showed the comparison between the predicted and the required cross sections for the unphysical pole 
for a given parameter set (Left). Here, we numerically solved the Boltzmann equation of the dark matter density. 
The figure shows that the predicted cross section satisfy the required value at J ~ 10"'' for Adyn = fn ~ 30 TeV. For 
larger values of Adyn and fa, the predicted cross section is always smaller than the required one. 

With the above constraint on Adyn and fn, we may derive an upper bound on the effective boost factor in the 
Breit-Wigner enhancement in the present model. The effective boost factor is given by, 

BP = ^, (49) 
and, for a give value of 7s, the factor is constraint from above by, 

BF < ^ ^ ^ . 30 X f ^) f-^V f . (50) 

Xf-K^s a;/r7i^ V^// \30TeV/ \ / 

As a result, we found that the effective boost factor cannot be significantly larger than 0(10) in this Nambu-Goldstone 
dark matter model. 

Second, we consider the case with the physical pole, S < 0. In this case, both the predicted and the required cross 
section decreases in in a region of |i5| > 7s. This means that the predicted cross section satisfies the required 
value only for a special value of Adyn and fa, i.e. 

50TeVx (^1 + ^^ . (51) 

Inversely, once Ajyn and fa take these values, the predicted dark matter density is consistent with observation for 
a wide range of S. In the right panel of Fig.[Tl we compared the predicted and the required cross sections. The 
figure shows that the predicted cross section is consistent with the required cross section in a wide range of S for 
Ady„~/fi~30TeV. 

Put it altogether, the both scenarios require the dynamical scale and the decay constant around Adyn — /fl — 30 TeV. 
From the analysis in the previous sectionllll these leads to the SUSY breaking scale F ~ A^/^ x 30 TeV. It should be 
noted, however, that there are 0(1) ambiguities associated with the strong dynamics to relate the dynamical scale and 
the SUSY breaking scale, and hence, the SUSY breaking scale can be slightly different from this naive expectation. 



See Ref . l26l for protohalo constraints on the effective boost factor of tlic dark matter annihilation. 
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TABLE I: Setup of the DM decay and annihilation 



Mode 




rria 


Decay of a 


Lifetime/Boost factor 


Decay 


5 TeV 


250 MeV 




1.5 X 10^" sec 


Decay 


1.5 TeV 


5 MeV 




5 X 10^" sec 


Annihilation 


2.5 TeV 


250 MeV 




1500 


Annihilation 


750 GeV 


5 MeV 




150 



III. COSMIC RAY SPECTRA IN NAMBU-GOLDSTONE DARK MATTER 

In this section, we discuss the cosmic ray spectra in the Nambu-Goldstone dark matter scenario. As we have seen, 
the right amount of the dark matter density is achieved by the annihilation process via the flaton resonance in the 
early universe. The cosmic ray spectra, on the other hand, depend on how the dark matters behave at the later time 
of the universe. In the Nambu-Goldstone dark matter scenario, we have two options for the behaviors at the later 
time. One option is the scenario with the enhanced annihilation and the other is the decaying dark matter. In the 
case of the annihilating dark matter, the annihilation cross section at the later universe is important. As we have 
seen, the unphysical flaton pole results in the Breit-Wigner enhancement, which greatly reduces the necessity of the 
astrophysical boost of the annihilation cross section. In the case of the decaying dark matter, the cosmic ray spectra 
can be explained regardless of the annihilation cross section, and hence, scenario works well for both the physical and 
unphysical pole scenarios. 

In the foUowings, we consider both options. As we will see, both options can fit the observed positron/electron 
spectrum. We further discuss the distinctive prediction on the gamma-ray spectrum, which comes from the finite 
annihilation/decay rates into the SSM particles. 

J^. The electron /•positron excesses from the Nambu-Goldstone dark matter 

We first consider the electron/positron spectra from the decay/annihilation of the Nambu-Goldstone dark matter 
with the final states consist of a pair of the R-axions which subsequently decay into pairs of electrons or muons.^^ 
The detailed parameter sets used for the decay/annihilation of the dark matter are given in TableU 

In Fig. [51 we show the predicted positron fraction (left) and the electron plus positron total flux (right). For 
the analysis on the propagation of the cosmic ray in the galaxy, we adopt the same set-up in Ref. [30j based on 
Refs. [Sll, [131, namely the MED diffusion model [SJ] and the NFW dark matter profile [s^l- As for the electron and 
positron background, we borrowed the estimation given in Refs. [H, [s^, with a normalization factor fcbg = 0.65. 
In the analysis of the positron fraction, we have taken into account the solar modulation effect in the current solar 
cycle [sgJ, which affects the fraction in £' < 10 GeV. The figures show that the excesses observed in the PAMELA and 
FERMI experiments are nicely fitted by both the annihilation/decay scenario with either four electron or four muon 
final states. 



5. Flaton decays into SSM particles 

In the above analysis, so far, we have considered that the dark matters mainly annihilate/decay into pairs of R- 
axions which eventually decay into muon or electron pairs. In the Nambu-Goldstone dark matter scenario, however, 
the interactions between the flaton and the SSM fields are not arbitrary but are fixed for a given gauge mediation 
mechanism. Thus, in the Nambu-Goldstone dark matter scenario, the SSM final states of the dark matter annihila- 
tion/decay are definitely determined for a given gauge mediation model, which gives the model a distinctive prediction 
on the gamma-ray spectrum. 

The decay modes of the flaton into the SSM particles are given as follows. Since we are assuming the models with 
gauge mediation throughout the paper, the interactions between the flaton and the SSM fields are obtained along 
with the gauge mediation effects. For example, the effective coupling between the flaton and the gauginos is given by 



Once the muon modes of the R-axion is open, the electron mode is ncgUgible, and hence, we do not need to consider the final state with 
the mixed leptons e+e~/^+/i~. 
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FIG. 2; Cosmic ray signals in the present model, (a): positron fraction with experimental data [ll. [27l. [2^ . (b): positron and 
electron fluxes with experimental data 0, 0, 0, ■ The yellow zone shows a systematic error and the dashed line shows the 
background flux. 



a Yukawa interaction; 



2 OS 



(52) 



where rrii denotes the gaugino mass and i runs the SSM gauge groups. For example, in a class of the so-called minimal 
gauge mediation [lB| , we obtain, 



drrii rrii 
ds fn 



In generic gauge mediation models, it is expected 



drrii 



7. 



R 



(53) 



(54) 



with an order one coefficient Cg}^ 



The coupling between the flaton and the sfermions and Higgs bosons are also obtained in a similar way; 



^oS — 



f 



ds 



X sff, 



(55) 



s=0 



and the model dependent coefficient dm^-fds is again given by. 



dm'j 



ds 



f 



(56) 



for the minimal gauge mediation model. 

With these interactions, the flaton decays into a pair of the SSM particles. For instance, the decay rate into a pair 
of the gluinos are given by 



F. 



^2 / \ 2 3 



47r \ms J f 



(57) 



In a class of model based on the messenger sector (37ll . the coefficient Cg can be vanishing for a special value of fji (see Fig. 4 in Ref. |38|). 
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Energy [GeV] Energy [GeV] 

(a) Annihilation (b) Decay 

FIG. 3: Gamma-ray signals in the present model with experimental data (0° < £ < 360°, 10° < |6| < 20°) ,45j. I and II 
represent the SSM spectrum where m^o — 100 GeV and 450 GeV, respectively. FSR+ICS represents the contributions of the 
FSR and ICS of the electron and positron cosmic ray. 



Thus, when the gluino mass is close to that of the flaton, the branching ratio of the flaton into the gluino pairs can 
be sizable, which gives a non-trivial contribution to the gamma-ray spectrum. ^"^ 

In the rest of this section, we demonstrate how the SSM modes affect the cosmic gamma-ray spectrum by concen- 
trating on the effects of the sizable gluino branching fraction and take the gluino branching ratio, 

Br, = ll^ = 8cl(^]\ 5.6 x 10^ x cl (^^Y ( > (58) 

as a free parameter, which we assume is typically below 10%.^"* The following analysis can be extended straightfor- 
wardly to more detailed analysis on model by model basis for a given gauge mediation mechanism, although we do 
not peruse in this study. 



6. Gamma-ray and antiproton spectrum with a gluino mode 

The gamma-ray signals come from the final state radiation (FSR), inverse Compton scattering (ICS) and the 
fragmentations of the final states of the decay/annihilation of the dark matter. Especially, in the present model, the 
contribution from the final states with gluinos has characteristic feature in high-energy region. (As for the FSR and 
ICS contribution, see Refs. HO, El and (S^. |43| . respectively) 

In this study, we concentrate on the gamma-ray signal from the fragmentation of the gluino final states. As a 
demonstration, we assume that the branching ratio to the gluinos of the dark matter decay/annihilation is 5%. 
Notice that the gamma-ray signal also depends on the decay mode of the gluino. Here, we assume that the sfermions 
are heavier than the gauginos, and consider two cases with m^a = 100 GeV and 450 GeV, while the gluino mass 
is fixed to mg = 500 GeV. For simplicity, we further assume the case that the gluino decays into only light quark 
pair q,q and the neutralino Xi^ the neutralino Xi decays into a gravitino and photon. In Fig. [3l the gamma- 
ray fluxes are shown. To estimate the fluxes, we have used the NFW proflle and averaged the halo signal over the 
region 0° < ^ < 360°, 10° < |6| < 20° For the annihilation case, we only include the halo component. In Fig. O 
we also show the FSR and ICS contributions from the electron and positron which come from the DM dominant 
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Notice that the branching ratios into the gauge boson pairs are highly suppressed |39| . 
In the decaying dark matter scenario, ms in E g. 1158 l l is replaced by m^. 
As for extra-galactic gamma-ray, see e.g., Ref. [44| 
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Energy [GeV] Energy [GeV] 



(a) Annihilation ttidm ~ 750 GeV. (b) Annihilation ttidm = 750 GeV. 




Energy [GeV] Energy [GeV] 

(c) Decay ttidm = 1500 GeV. (d) Decay moM = 5000 GeV. 

FIG. 4: Antiproton signals in the present model with experimental data [4^. I and II represent the SSM spectrum where 
m^o = 100 GeV and 450 GeV, respectively. 



annihilation/decay mode. We estimate the ICS component with the method discussed in Ref. [43i| . using data of 
interstellar radiation field provided by the GALPROP collaboration [4^ , which is based on Ref. |47l|. 

In both cases, the DM signals are consistent with the current experiment data, and anomalous behavior of the 
gamma ray is expected around the DM mass for the annihilation cases or half for decay. This behavior comes from 
the hard component from the neutralino Xi decay. 

Th e g luino final states also raise the antiproton signal, which is severely constrained by the PAMELA experi- 
ment [45. In Fig. m we show the ratio of antiproton which comes from signal and the background proton. In this 
analysis, we used the program PYTHIA 49] to obtain the fragmentation function of the final states into antiprotons. 
The analysis of the propagation of the antiproton is again based on that used in Ref. [30j] . The figure shows that the 
antiproton fiux strongly depends on the diffusion models, and the shaded region corresponds to the dependence on 
the diffusion model, with upper side of the region corresponds to the diffusion model MAX, under side MIN and the 
line MED in Ref. [33] . The proton background is taken from Ref. [50] . The larger antiproton fluxes in the case of "I" 
reflect the higher energy fractions of the quark pairs in the decay of the gluinos. We see the contradiction between 
the experiments and the signals in some diffusion models even for the 5 % branching ratio into the gluino final states. 

As a result, we see that the antiproton flux provides very strong constraint on the model. For example, in the 
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minimal gauge mediation, the branching fraction into the gluino final states are solely determined by the masses of 
the dark matter and the gluino, i.e. Cg = 1 in Eq. (j58p . Thus, the strict constraint on the branching ratio leads to a 
strict constraint on the mass of the gluino which is an important parameter for the SUSY search at the coming LHC 
experiments. Therefore, the Nambu-Goldstone dark matter scenario can be investigated with the interplay between 
the cosmic-ray experiments as well as the direct SUSY search at the LHC experiment. 
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